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ABSTRACT 

Aims. Binary interactions are believed to be important contributors to the structures seen in planetary nebulae (PN), and the sole cause 
of the newly discovered compact dust disks. The evolution of such disks is not clear, nor are the binary parameters required for their 
creation. 

Methods. We study the Galactic bulge planetary nebula M 2-29 (for which a 3-year eclipse event of the central star has been attributed 
to a dust disk) using HST imaging and VLT spectroscopy, both long-slit and integral field. 

Results. The central PN cavity of M 2-29 is filled with a decreasing, slow wind. An inner high density core is detected, with radius 
less than 250 AU, interpreted as a rotating gas/dust disk with a bipolar disk wind. The evaporating disk is argued to be the source of 
the slow wind. The central star is a source of a very fast wind (~ 10' kms"'). An outer, partial ring is seen in the equatorial plane, 
expanding at 12kms"'. The azimuthal asymmetry is attributed to mass-loss modulation by an eccentric binary. M2-29 presents a 
crucial point in disk evolution, where ionization causes the gas to be lost, leaving a low-mass dust disk behind. 

Key words. ISM: planetary nebulae: individual: M 2-29 (PNG 004.0-03.0) - Stars: AGB and post-AGB - planetary nebulae: general 



1. Introduction 

Disks are a recent addition to the range of structures seen in plan- 
etary nebulae. The outer nebulae show a large variety of morpho- 
logical characteris tics, including tori, bipola r flows, and ellip- 
soidal shapes (e.g. lRamos-Larios et al.ll2008h . with typical sizes 
of lO^'-lO^ AU. But embedded in these nebulae, in a number of 
cases a central dust disk has been discovered with sizes of or- 
der 1 0^ AU, which can be resolved with the VLT interferometer 
(e.g. IChesneau et al ] l2007h . The disks appear well aligned with 
the outer nebulae but are too low mass to directly have affected 
their formation. More likely they are a byproduct of the original 
shaping process. 

Other evidence exists for the presence of material very 
close to s ome central stars, inc luding unresolved emission- 
line core s (Rodri guez et al.l l2001) and eclipses from dust disks 
dHaiduk et al. 2001). The disks are believed to require binary 
systems for their formation but the details of the required system 
parameters are unclear. An interesting review on binary syst ems 
and th eir shaping of PNe was recently published by de Marcol 
(l2009l) . 

Here we present a study of a planetary nebula with such an 
unresolved emission-line core, known to harbour a dust disk: 
M2-29 (PNG 004.0-03.0). Deep echefle spectra and integral- 
field spectta, together with HST imaging, allow us to disentangle 
the various components of this complex nebula, and to determine 
the velocity fields. We find the central core to be the source of 



* Based on observations collected at the European Organisation for 
Astronomical Research in the Southern Hemisphere, Chile (proposals 
075.D-0104, 077.D-0652, 081.D-0130) and HST (program 9356) 



an on-going, but decreasing wind, which we attribute to a disk 
wind. Although direct evidence for the binary companion re- 
mains elusive, the outer structure suggests an interaction with a 
star on a non-circular orbit. From comparison with other objects 
with dust disks, we find evidence that such disks strongly de- 
crease in mass during the planetary nebula evolution. However, 
they remain sufficiently massive that residual disks can survive 
long into the white dwarf phase. 



2. Observations 

2.1. The HST images 

SNAPshot HST/WFPC2 images of M 2-29 were obtained in 
2003, in three diff^erent filters: F656N, F547M and F502N. 
The pixel scale of the camera is 0.0455 arcsec. The Ha image 
(iHaiduk et al.ll2008h is show in Fig.[T] The main nebular compo- 
nents are the central source, a wing-hke structure, and the ex- 
tended, fainter nebula. 

2.2. VLT/UVES long-slit spectra 

A 600-second VLT/UVES (iDekker et al.l2000l) echefle spectrum 
was obtained in 2005, with the slit 0.5 arcsec wide and 1 1 arcsec 
long. The slit was put through the centre of the nebula, aligned 
with the jet-like structure. The pipeline calibration was applied, 
including wavelength calibration and merging of the separate or- 
ders of the spectrum. The observations cover a spectral region 
approximately from 3300 to 6600 A. The long-slit spectra of 
M 2-29 are spatially resolved although they have much worse 
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resolution than HST images and the seeing during the observa- 
tions was nearly 1.5 arcsec. 

2.3. VLT/ARGUS-IFU spatially resolved spectroscopy 

We also retrieved publicly available data from the ARGUS 
Integral Fie ld Unit (IF U), a mode of the FLAMES instrument 
at the VLT dPasquini et al. 2002). These observations of M 2-29 
were performed in visitor mode as a part of program 077.D-0652 
(PI. D. Schonberner). Two spectral ranges are available, HR8 
and HR14B covering range 4911-5158 and 6383-6626 A with 
resolution R~32000 and R~46000 respectively. Because we 
aimed for small scale structure, we selected only those frames 
which had seeing better than 1 arcsec; this constraint yielded 
one frame in HR14 (700s) and six frames in HR8 (6x1 10s). The 
large-area (12x7 arcsec^) IFU was used, where each single re- 
solving element (spaxel) covers 0.52x0.52 arcsec-. 

The data were analyzed with the girBLDRS pipeline (a.k.a. 
Geneva pipeline). Correction of bias level flat fielding and the 
wavelength calibration were prepared in a standard way. Spectra 
were extracted with the sum-extraction method, because the de- 
fault optimum extraction is not well suited to emission-like ob- 
ject. The airmass was lower than 1.06 for all frames alleviating 
the need to correct for differential atmospheric refraction. The 
data were not flux calibrated. 



2.4. VLT/VISIR images 

M 2-29 is a known dust emitter, detected at 25 jum by IRAS (the 
12-/im detection was rejected for the faint-source catalogue, per- 
haps because of con fusion with a nearby star). N-band spec- 
troscopy presented bv lCasassus et alj ( 1200 lb indicates a possible 
broad silicate emission band at 8-12//m. We retrieved and re- 
duced archival Spitzer IRS data, which confirms the presence of 
silicate emission, and also shows a weak PAH feature at 1 1 .2 fim 
(Fig.H. 



M2-29 
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. 1 ■ 









Fig. 2. Left: Spitzer low-resolution spectrum, showing a broad 
silicate band with superposed emission lines and a weak 11.2 jum 
PAH band. Right: the VISIR image showing a strong point 
source. The image is 3.8 x 3.8 arcsec in size. The transmission 
curve of the filter used is shown on the spectrum. 



To study the location of the dust, we obtained mid-infrared 
images of M 2-29 using VISIR on the VLT (Lagage et al. 2004). 
We used a 0.075 arcsec per pixel scale and and the VISIR SiC fil- 
ter (Ac - 11.85/im, AA = 2.34jum). The observations were car- 
ried out in visitor mode with great weather conditions, leading 
to diffraction limited images during the whole run. The observa- 
tions were done using the standard chopping and nodding tech- 
nique to reduce the background emission in the mid-infrared. 
The data reductio n was performed using self-developed IDL rou- 
tines described bv lLagadec et al ] (l2008h . Images were corrected 
for bad pixels and then co-added to produce a single flat-field- 
coiTected image, comprising the average of the chop and nod 
differences. 

The resulting image at 1 1 .85 fim shows an unresolved source 
(Fig.|2]i, indicating that the emitting source is smaller than ~0.3 
arcsec. Thus, the emitting hot dust only comes from one compact 
component (the core), and not from the extended nebula or wing- 
like feature. 



K. Gesicki et al.: Disk evaporation in a planetary nebula 



3 



-a 

0) 

o 
m 

'm 

1) 



Ha cross sections 

along wing 

_ _ perpendicular 




1 

0.8 
0.6 
0.4 
0.2 




Fig. 3. The cross-sections through the H a HST image of M 2- 
29 taken along the wing-like structure and perpendicularly to it. 
The intensity scale is the same for both profiles 



3. Morphology from HST images 



The HST Ha image dHaiduk et al.1 |2008|) shown in Fig.[T] 
presents the complexity of the nebula, including the main ring, 
the halo, the wing-like structure, and the unresolved central ob- 
ject. The main ring is actually composed of two arcs, which may 
not quite connect. The gaps between them are nearly perpendic- 
ular to the wing-like structure. There is a slight indication for 
elongation or brightening in the EW direction, beyond the gaps 
(best seen in the left panel, with a linear scale). The overall im- 
pression is that of a barrel-shaped nebula, with a beginning bipo- 
lar flow, a relatively common shape among planetary nebulae. 

The HST Ha and [O iii] 5007 A images are almost identical 
and show the same extent. This means that M 2-29 should be 
fully ionized (density bounded). For the discussion below we 
use the Ha image. 

3.1. The outer halo and dense wind 

To inspect the outer regions, we overplot in Fig. [3] two scaled 
cross-sections: one extracted along the wing-like structure and 
the other perpendicularly to it. The outermost nebular regions 
decrease in brightness in the same way in all directions, indica- 
tive of sphe rical symmetry. The re is no clear evidence for ISM 
mteraction dWareing et al.i2007h . 

Intriguingly, the emission profile extends smoothly and sym- 
metrically inside the arcs, up to the central object. This is shown 
in (Fig.O, where the small bump on the left hand side corre- 
sponds to the arcs. The arcs and jet appear to be superposed on 
a largely spherical component whose brightness smoothly de- 
creases outward. This component we will call the "dense wind". 
The spherical symmetry of this component points to a stellar 
origin. Its brightness at Ha image proves it is not hydrogen defi- 
cient. 

The emission profile of the dense wind, including the halo, 
corresponds roughly to the surface brightness of a sphere filled 
with matter of a r""'^ density distribution. This is unusual. The 
m ost widely assumed A GB wind density is r^^, as discussed 
by IVillaver et alj ( l2002h . The r density implies a mass-loss 
rate decreasing as M oc r^-^. The surface brightness with an ex- 
pansion velocity of lOkms"' gives a current mass loss rate of 
lO^'^Moyr"'. 
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Fig. 4. The cross-sections through the Ha HST image of M 2-29 
taken along the wing-like structure. The dashed line shows the 
observed intensity profile, the dotted lines show both spherical 
components - the arcs (photoionization model) and the dense 
wind, the continuous line presents the difference between the 
observed cross section and the two spherical components. 



3.2. The central nebula 



From the HST images, dHaiduk et al.ll2008l) showed that the un- 
resolved source at the centre contains an ionized nebula which 
outshines the central star in the emission lines. 

A number of artificial P SFs were created using the Tiny Tim 
v.6.3 dKrist & Hookl l2004 . and fitted to this central nebula. 
Series of sub-pixel images for different shifts relative to the 
nebular image were obtained. The PSFs were convolved with 
Gaussian kernels of different FWHMs, and re-sampled to the 
resolution of the Planetary Camera. After re-sampling, the 
charge diffusion kernel was applied. The PSF was then sub- 
tracted from the M2-29 image. Minimum residuals were ob- 
tained for PSFs convolved with a Gaussian of FWHM ~ 0.8 
pixel. A similar result was obtained for a background star The 
broadening of the PSF with respect to the theoretical one could 
be due to an imperfect focus. From the comparison with the field 
star, we constrain the deconvolved FWHM of the nebula as < 0.6 
pixel. This coiTesponds, at the distance of 8 kpc, to a diameter of 
< 250AU(3x 10'^ cm). 



3.3. The wing-like feature 

To obtain the radial emission profile of the wing-like feature 
from the HST image, the other nebular components should be 
subtracted. 

The profile of the dense wind was obtained from the dashed 
line in Fig.|3] by taking the average of this cross section with 
its mirror reflection. The arcs are more difficult, as they are not 
circular symmetric and cannot be taken from the perpendicular 
profiles. Instead we used a photo-ionization model (see Section 
5), approximated by a spherical shell best fitted to the two arcs, 
and subtracted the profile calculated from this model. The unre- 
solved core was not subtracted. At the decomposition weighting 
factors were applied. The result is shown in Fig.|4] 

To the left (south) of the central object, a flat profile was ob- 
tained. The wing-like structure shows up only to the north (right) 
of the core. 
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No clear residual is seen at the cross-over point of the wing- 
like feature and the arcs. A superposition on the sky, rather than 
a physical connection, seems to be indicated. 

3.4. The arcs 

The successful decomposition of the image cross-section into 
the dense wind and the photo-ionization model, shown in Fig.|4] 
indicates that the southern arc can be interpreted as a fragment 
of a photoionized sphere showing emissivity enhancement at the 
edge. The same can be said about the northern arc. A very com- 
mon interpretation of such a structure is that the arcs form the 
projection of a thin barrel, open to the two poles. The dense wind 
fills the inside and the outside of the barrel. This appears to sug- 
gest that the density enhancement of the barrel formed out of the 
dense wind. 



4. The spectra 

4.1. Velocity maps 

In Fig.|5]we present spectroscopic observations in emission lines 
of [Om]5007A and [Nn]6583A. These are strong forbidden 
lines, which probe the inner and outer nebular regions respec- 
tively. For both lines we extracted from the data cubes three rep- 
resentative slices corresponding to the systemic radial velocity 
(panels b and ffl and those shifted by 12kms ' below (panels 
a and e) and above (panels c and g) this value. The grey scale 
represents the relative intensities normalized for each slice sepa- 
rately. The square pixels are of spatial size of 0.52 x 0.52 arcsec^. 
The contours show the HST images; [Om] for the upper row. 
Ho- for the bottom row. The right-most frames show the long-slit 
UVES spectra, where the slit was positioned along the wing-like 
structure i.e. 5° away from the N-S direction. 

At the central star position, the [O m] velocity maps show 
strong emission at both the positive and negative velocities (pan- 
els a and c) but a minimum at the systemic velocity (panel b). 
This emission corresponds to the unresolved central component 
and can be interpreted in terms of a bipolar outflow. The wing- 
like structure can be traced in panels (a) and (b) with velocity 
decreasing outwards; this can also be seen in the long-slit spec- 
tra (panel d). The southern arc is best distinguished at the zero 
velocity panel (b). 

More details can be derived from [N ii] images. The wing- 
like structure at -12kms"' (panel e) is clearly more extended 
than in [O iii] at the same velocity (panel a), and is seen at both 
sides of the central component. The systemic velocity image 
(panel f) shows the wing-like structure extending two arcsec- 
onds in the northern direction, far beyond what is seen in the Ho- 
image. The long-slit UVES spectra (panel h) also show this sig- 
nificant extension. The southern arc can also been at the systemic 
velocity, and a very weak component is seen in panel (f) extend- 
ing for a couple of pixels in S-W direction. Panel (g) shows that 
a component at the eastern end of the southern arc is moving 
away from us, in fact at a little more than -i-12kms 

In Fig|6]we present more examples of the 2-D (position- 
velocity) contours presentation of bright emission lines in the 
UVES spectrum. All plots are shown to the same scale, but in the 
blue spectral region the slit was shorter so the spatial extent of 



' The veloci ty is corrected for the radial velocity of the object of 
-11 2.2 km s"' dPurand et all 199^ and barycentre velocity for the mo- 
ment of observations of - 10.4 km s"' for UVES and -22kms"' for 
ARGUS. 
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Fig. 7. Model intensity distribution for Hi, Om and Nil The 
overall density follows H to a good accuracy. 

Table 1. Dereddened M2-29 emission line flux ratios, relative 
to I(Hy6)=100.0. Observed d ata come from our S A AO spectrum 
(dereddened with C/? = 1 . 1 8 (Hai duk et aDl2008h ) and from liter- 
ature (dereddened by the authors). 



Wavel. 


Ident. 


Observed Flux 




Model 






SAAO 


(1) 


(2) 




3726 


[On] 


43.6 


42.1 


34 


15.8 


3869 


[Ne in] 


72.2 


49.8 


58 


19.7 


4363 


[Oiii] 


27.0 


14.0 




1.9 


4959 


[Oiii] 


150.2 


146.0 




46.8 


5007 


[Om] 


449.4 


424.0 


452 


134.9 


6300 


[Oi] 


1.6 


1.1 




0.4 


6548 


[Nil] 


7.0 


6.9 




10.2 


6584 


[Nil] 


17.8 


20.8 


17 


29.9 


6716 


[Sii] 


1.2 


1.9 




0.7 


6731 


[Sii] 


1.8 


2.4 


2 


0.8 


7136 


[Ar m] 


9.0 


12.5 




8.3 



References. (l) Exteretal. (2004) ; (2) Howard et al. (1997) . 



these lines is smaller. Many more emission lines are identified in 
the spectrum, but the three presented here best supplement those 
in Fig.|5] showing the variety of line profiles. 

The emission of the wing-like feature is visible from the cen- 
tre until about 2 arcsec up. It is very pronounced in the [O iii] 
5007 A, [Nil] (Fig.|5]l and Hei lines, very weak in the [Oi] line 
and completely absent in the highly excited [O iii] 4363 A line. 

The velocities are relatively low (around 10 kms"') com- 
pared to typical expansion velocities of planetary nebulae. 

4.2. Photoionization modelling of the outer nebu/a 

To model the spectra a nd line profiles, we apphed the Torun 
photo -ionization codes dGesicki & Ziilstral 120031: iGesicki et al.l 
I2006h . The star is assumed to be a black body with a luminos- 
ity and effective temperature. The nebula is approximated as a 
spherical shell defined by the radial density distribution and the 
radial velocity field. The chemical composition ([ O/H] =7.44) 
and th e observed line intensities were adopted from lExter et al.l 
(l200l . We have also measured the line intensities from the 
low resolution SAAO spectrum dHaiduk et alJ l2008h . TabIe[I] 
presents the observed data together with our best model. The 
model emissivities are shown in Fig|7] 



K. Gesicki et al.: Disk evaporation in a planetary nebula 



5 






(a) [Oin]5007A at -12kms 



(b) [Oiii]5007A at system velocity (c) [Oiii]5007A at +12kms 







(d) [O III] 5007A long-slit spectrum 




(e) [Nii]ti583A at-12kms 



(f) [Nii]6.583A at system velocity (g) [Nii]6583A at +12kms 



-40 -20 20 40 
uelacity [km/s] 

(h) [N ii] t)5S3A long-slit spectrum 



Fig. 5. The ARGUS-IFU and UVES ob- 
servations of M2-29 in spectral lines of 
[O in] 5007 A (upper row) and [Nn]6583A 
(lower row). The presented slices (left) are 
extracted from ARGUS data cubes at wave- 
lengths corresponding to velocity shifts from 
the systemic velocity of -12kms"', and 
-1-12 km s"' (from left to right respectively). The 
fluxes are shown in grey in linear scale. The 
black colour corresponds to the relative inten- 
sity of 0.9 except of panels (e) and (f) where 
cuts were applied at levels 0.5 and 0.3 respec- 
tively. In both lines the image at -12km s"' 
(panels a and e) is the brightest. For [O in] the 
central image (b) is of similar brightness while 
(c) is two times weaker. For [N n], panel (f) and 
(g) are fainter by a factor of 3 and 30 with re- 
spect to panel (e). The contours are from the 
HST images, overplotted to the same spatial 
scale; for the [Oin] slices the [Om] image is 
shown while for the [Nii] slices the Ha im- 
age is shown. The right-most panels show the 
UVES long-slit spectra where the slit was posi- 
tioned along the wing-like structure i.e. nearly 
N-S direction. The contours are 0.004, 0.02, 
0.1, 0.5, 0.95 for [On] in panel (d) and 0.007, 
0.03, 0.08, 0.3, 0.8 for [Nn] in panel (h). 




3 -20 20 40 
ijelocity [Km/s] 
(a) Hei5876A 




-4 - 



-40 -20 2C 40 
velocity [km/s] 



(b) Oi6300A 



10 -20 20 40 
^^e]QCity [km/s] 
(c) [Ora]4363A 



Fig. 6. The VLT/UVES long slit spectrum of 
He I 5876 A (left), [Oi] 6300 A (middle) and 
[O m] 4363A (right). The velocity scale is cen- 
tred on the systemic velocity. The wing-like 
structure extends up from the centre. The con- 
tours are at levels of 0.02, 0.07, 0.2, 0.7, 0.9 in 
panel (a), 0.06, 0.1, 0.2, 0.5, 0.9 in panel (b) 
and 0.01, 0.03, 0.1, 0.3, 0.8 in panel (c). 



The one-dimensional model by necessity assumes spherical 
symmetry. It solves for the observed line intensities and the ra- 
dial surface brightness distribution. For the latter, we extracted 
the image cross-sections taken at 45 degrees left and right to 
the wing-like structure i.e. avoiding the wing while including 
the arcs. The two cross-sections are quite similar The model as- 
sumes a central cavity: neither the central unresolved nebula, nor 
the inner part of the dense wind (i.e. inside the arcs) is included. 

The satisfactory solution is summarized in Table|2] Although 
the comparison of the Ha and [Oiii] HST images indicates a 
density bounded PN, the line intensities (integrated over the 
whole nebula) are in fact significantly better fitted when the PN 
is ionization bounded. We solved this problem by adopting the 
density distribution extending far beyond the arcs, obtaining the 
ionization front at about 4 arc sec radius. Th i s mak es the model 
slightly different from that in iHaiduk et al.l (l2008h . The strong 
[O m] 5007 A line is not well fitted. We could only fit this line 
by increasing the oxygen abundance by a factor of 2, still metal- 



poor however we did not attempt to improve the abundance de- 
termination. 

The model shows a constant expansion velocity with dis- 
tance from the star, as derived from the similar widths of lines 
formed at different radii, broadened by turbulence. This is at 
variance with the usual positive velocity gradient seen in plan- 
etary nebulae. The model turbulence probably reflects the su- 
perposition of different nebular components. The mass-averaged 
expansion velocity derived for M 2-2 9 is 12kms~' (cornpared to 
14kms ' obtained from [Oiii] only. iGesicki & Ziilstral hoOOh ). 
Assuming a distance of 8 kpc this velocity results in a kinematic 
age of about 5000 yr 

The ioniz ed mass of the nebula i s rather high, in view of the 
conclusion of iTorres-Peimbert et a concerning old age 

(and therefore low initial mass) of the object. The age may have 
been overestimated, or mass transfer can be proposed. However, 
caution should be applied when deriving a mass from a spherical 
model for a non-spherical nebula. 



6 



K. Gesicki et al.: Disk evaporation in a planetary nebula 



Table 2. M 2-29 photoionization and kinematic model parame- 
ters. 



Parameter 


Value 


distance 


8kpc 


^cff 


70 000 K 


Luminosity 


1200 L© 


Ionized Mass 


0.38 Mq 


log He / H + 12 


10.99 


log N / H + 12 


7.14 


logO/H+ 12 


7.44* 


logNe/H + 12 


6.84 


log S / H + 12 


5.64 


Nebula Expansion Velocity 


12kms-i 


Turbulent Velocity Component 


7 km s"' 


Kinematic Age 


5 000yr 


* this adopted value might be too low (see the text) 




-600 -400 -200 200 400 600 



1,^6 1 DC i ty [ Km/B ] 

Fig. 8. The P-V plot of the Ha line observed with UVES, with 
the long-slit aligned with the wing-like structure. At the right 
edge the permitted line of C ii 6578 A can be seen. Contours are 
equally spaced in logarithmic intensity. 

4.3. Extended wings of Ha 

The Ha emission line reveals very wide wings, seen at the 
location of the central unresolved source only. In the UVES 
position-velocity diagram (Fig. [8]), the vertical structure shows 
the emission from the wing-like structure, and the horizontal ex- 
tent shows the high velocity wings, located at the same position 
as the stellar continuum. 

Averaging the slit spectra over the central 1.5 arcsec, the 
emission wings can be traced over 1000km s"' (Fig.|9]l. The 
ARGUS-IFU data cube confirms that at wavelengths beyond 
+4 A ( +200 km s"' ) from the line centre the emission originates 
solely from the unresolved source. 

Extended wings can also be seen in [O iii] 5007 A shown 
in Fig.|5}i), and in Hfi (not shown here). For weaker lines it is 
not clear whether the wings are present or absent (see e.g. N ii 
6583 A line in Fig.|5h). However, the Oiii5007A line shows 
wings one order of magnitude narrower than Ha. This should 
be considered when interpreting this emission component. 

The fact that the large broadening is only seen in some lines 
argues against electron scattering. The presence of a high veloc- 
ity component in the central object is indicated. 




J \ \ \ I \ \ \ \ \ \ I \ L 

-1000 1000 

velocity [km/s] 



Fig. 9. The Ha line extracted from the central pixels of UVES 
(solid line) and ARGUS-IFU (dotted line). The emissions are 
normalized to the line peak. Note that ARGUS has a lower spec- 
tral resolution. 




Fig. 10. The ARGUS-IFU collapsed images of the Cii6578 A 
and He 1 5016 A lines. North is on top, east to the left, the size of 
a single pixel is 0.52 arcsec. 



4.4. The permitted emission lines 

The UVES spectrum shows a number of permitted emis- 
sion lines. The recombination lines of Cii at 4267.26 and 
6578.05 A are present (e.g. Fig.|8]l. Also the recombination lines 
of Om at 4638.85, 4641.81, 4649.14 and 4650.84A can be 
identified. These l i nes ar e useful in abundance analysis, e.g. 
iGarcia-Roias et al.l ( l2005h . 

For the well exposed C n 6578 A line, an ARGUS-IFU im- 
age was obtained by integrating over its line width. This line 
originates in an extended region as can be seen in Fig.fTOlwhere 
it is compared with He 1 5016 A line. Concerning the other per- 
mitted lines the situation is not so clear since they are weaker 
and more noisy. The well-known group around 4650 A falls at 
an echelle order merging where distortions may occur; it is also 
not covered by the ARGUS-IFU spectral range. The UVES data 
suggests that it may also form in an extended region but this re- 
quires confirmation. 

We also identified a number of permitted N iii lines (4097.33, 
4103.43, 4634.14, 4640.64, 4641.85 A). They ar e often inter- 
preted in terms of the Bowen mechanism (see e.g. lSelvelli et al.l 
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Fig. 11. The interstellar absorption line profiles of ions Na i and 
CaiiK in the spectrum of M 2-29. The profiles are shown in the 
heliocentric reference frame. 

(I2007h l where He ii Lya overlaps with an O m line, and the re- 
sulting O III cascade itself has a frequency coincide with an N iii 
resonance line. The O iii Bowen lines below 4000 A fall within a 
very noisy area of our VLT spectrum and can not be identified. 
The very weak He ii l ine argues again st a Bowen mechanism in 
the nebula of M 2-29 . iMihalasI d 1 97 ih proposed dielectronic re- 
combination for explaining the N m line complex near 4650 A. 
The N III lines are only detected from the central nebula. 

4.5. Interstellar absorption and distance 

The UVES spectra of M 2-29 show interstellar absorption lines 
of Na I and Ca ii, shown in Fig.fTTI 

The equivalent width (EW) of Ca ii H and K lines (i.e. adding 
together all Doppler components) shows an approximate linear 
relation with distance (Megier et al. 2005), as derived from OB 
stars with Hipparcos parallaxes. For M 2-29, the measured EW 
of the full CaiiK absorbing complex is about 1500 mA. (The 
H-band is confused with He.) This result s in a distan c e of 4 
or 6 kpc when applying the relations of iMegier et al.l ( l2005h 
orjOalazutdinov (2005) respectively. It is not clear whether the 
linear relation can be extrapolated into the inner regions of the 
Galaxy; however, the indicative values are in agreement with a 
Galactic Bulge location. 

The interstellar absorption line profiles show components 
from different interstellar clouds. The three main compo- 
nents are at (heliocentric) velocities of 0-50 km s"', at -60 to 
-20km s"', and at -100 to -80kms"'. The absorption at posi- 
tive velocities can be associated with the Molecular Ring around 
the inner Galaxy, while the intermediate negative velocity com- 
ponent is consisten t with the 3 kpc expanding ring as seen in CO 
(lDameetalJl2001h . However, the Galactic CO maps show no 
counterpart to the most negative velocity component in M 2-29. 
This component cannot be explained with Galactic rotation, nor 
Bulge rotatior0. 

HIPASS spectra in this region do show an Hi emission fea- 
ture corresponding to this negative velocity component. It may 
be related to a nearby catalogued high velocity cloud, HVC 
002.9-06.2-088 (Putman et al. 2002), with the same velocity, as 

- The mean rotatio n of the Galactic bulge at the longitude of M 2-29 
is around +50kms ' jMinniti & Zoccalill2008h . 



the emission at this velocity seems to extend over some dis- 
tance. The distance to this cloud or complex is not known, but as 
M 2-29 is clearly located behind it, we suggest it is situated rela- 
tively close to the Galactic centre. Small Ca ii and Na i absorption 
clouds with si milar high velocities ap pear to be common in the 
Galactic halo ( iBen Bekhti et al.|[2008h . with a variety of origins. 

It is remarkable that the non-rotational interstellar absorp- 
tion component near -100km s"' is very close to the systemic 
radial velocity of M 2-29. As a PN, or related object, M2-29 is 
expected to be old and evolved, and any association with the in- 
terstellar component is deemed accidental. 

5. Disentangling the nebula 

The nebula of M 2-29 shows four distinct components: the inner 
unresolved nebula, the dense wind extending outwards, the arcs, 
and the wing-like structure. The arcs appear to be an enhance- 
ment within the dense wind rather than an independent struc- 
ture. Here we attempt to obtain an understanding of the other 
two structures. 

5.1. An outer partial ring: the wing-like structure 

The 'wing-like' feature is a low-velocity structure. Based on the 
observed location and the lack of interaction with the arc, we 
locate it outside of the main nebula. We identify it as a partial, 
expanding ring, surrounding the main nebula, and whose plane 
is almost coincident with the line of sight. The ring has an inner 
radius of about 2 arcsec. 

The velocity structure along this ring is depicted in Figs. 
|5] and |6] In the [O iii] line, only the top (north) of the ring is 
seen, but the [Nii] line shows the more complete structure. In 
the southern part, the [Nii] is displaced to the left (see Fig.|5^). 
The displacement suggests that the plane of the ring is seen at an 
inclination of approximately 20 degrees with respect to the line 
of sight. Because the displacement is seen in the rear part of the 
ring it is not unlikely that the structure is warped. 

The line splitting shows that it is expanding: an expansion 
velocity (with a minor correction for the inclination) of about 
12 + 2kms"' is indicated. 

The UVES slit was aligned with the brightest part of this 
structure, however because of the inclination the southern ex- 
tent fell outside the slit. This is made clear by the ARGUS im- 
age. The position-velocity diagrams from the UVES data there- 
fore have an incomplete coverage. Still, the brightness distribu- 
tion of the ring is strongly asymmetric, while the velocities ap- 
pear symmetric. The peak [O iii] emission occurs at a velocity 
of -lOkms"', offset by about 1 arcsec from the star. This cor- 
responds to a location on the ring about 30 degrees from the 
line of sight to the centre. The bright [Oiii] emission occurs over 
an arc of roughly 90 degrees, on the forward side of the ring. 
Some faint emission originates from the backside, seen towards 
the southwest at positive velocities, in Fig.|5};- 

The bright [N ii] emission covers a much wider arc, extend- 
ing for some 120 degrees. At a much fainter level the ring may 
be almost closed. The [N ii] emission is spatially also much more 
extended than [Oiii], with a fainter outer region with a radius of 
at least 4 arcsec. This is seen in Fig.|5f as the northern extension. 
This region only contributes at near-systemic velocity, suggest- 
ing that the expansion velocity decreases sharply from 2 to 4 
arcsec radius. The radial extension is not seen towards the south, 
showing that at its outermost radius the ring is is not closed, per- 
haps limited to some 180 degrees. 
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The ring is also seen in Hei, C ii, and [Oi], at the location of 
the brightest [O m] and [N n] emission. Thus, there is sufficient 
column density in this direction to cause ionization stratifica- 
tion and an ionization boundary, as also indicated by the photo- 
ionization model. The lack of [Om] beyond 2 arcsec is consis- 
tent with the photoionization model of Fig.|7]. 

Finally, this ring system is located close to the equatorial 
plane of the main nebular barrel defined by the northern and 
southern arc (see Sect.|3]l. A slight misalignment is explained 
by the inclination. 

Thus, we arrive at a model of an expanding equatorial, par- 
tial ring, 2-4 arcsec in radius, with a strong azimuthal density 
gradient, a radial density gradient, and an expansion velocity of 
12 + 2kms"' at 2 arcsec, decreasing with increasing radius. 

5.2. The unresolved disk 

The central source is detected in emission lines of Ha, [O m] 
4363 and 5007 A, Hen, the Bowen Niii Unes, and [Oi]. Other 
lines are present but are confused by the spatial and veloc- 
ity overlap with the extended emission. The very high velocity 
wings, and the [O m] 4363 A line are only seen in this core. The 
size of the core is < 0.03 arcsec, or < 250AU at the distance of 
the Galactic Bulge. 

The [O III] 4363/5007 line ratio is much higher in the core 
than in the extended nebula0 T his indicates high dens i ty and 
temperature, in agreement with iTorres-Peimbert et al.l (Il997h 
who first showed the presence of a high density and high elec- 
tron temperature component close to the star. We ran some test 
models to estimate the parameters of the inner nebula, using the 
central star parameters derived above. In these models, the ratio 
4363/5007 [Om] lines can be reproduced with Te ~ 1.4 x 10"* K 
and We ~ 6 x 10^ cm~^. These values are a little lower than the es- 
timates of lTorres-Peimbertet a n (ll997h . At these densities, the 
mass of the compact nebula needed to reproduce the observed 
Hj8 flux is 10"^ Mq. This value is not very accurate, as the core 
is unlikely to show constant density and spherical symmetry, 
but provides a believable order of magnitude. Other lines often 
used as the electron density indicators, e.g. lower ionized/excited 
[O ii] or [S ii], predominantly originate in the wing-like structure. 
This makes an extraction of the central component much more 
unreliable, in no case as clear as the [O iii] 4363 A line. 

All lines show velocity broadening of ~ lOkms"'. For this 
velocity, the kinematical age of the core is less than 100 years, 
requiring a continuing mass-loss rate of ~ lO"^M0yr"'. The 
lower limit to the mass-loss rate is 1 to 2 orders of magnitude 
lower than that required for the extended nebula (assuming a 
kinematic age of 5000 yr), and the surface brightness distribu- 
tion of the extended wind indeed implies a linearly decreasing 
mass loss rate with time. It raises the question, however, why the 
mass-losing star is not seen. Such mass-loss rates and velocities 
are found in red giants. (The wind from the hot central star will 
be very much faster). Neither the spectrum nor the photometry 
shows evidence for a contribution from any other star than the 
PN central star The 2MASS photometry excludes the presence 
of a red giant more luminous than 100 L©. Although a mass- 
losing star has clearly been present in the recent past, it does not 
appear to be present at the moment. 

The alternative to on-going stellar mass loss is a stable, 
rotating structure. This has some evidence in favour: the line 



In the photoionization model of the outer nebula the [O m] 4363 A 
line is underestimated by an order of magnitude which is another argu- 
ment that it does not belong to the outer main ring. 



widths decrease with decreasing excitation, from +18kms"' 
HWHM for [Oiii] line to close to a few kms ' or less for the 
[O i] line consistent with rotation velocities decreasing outward. 
Hajduk et al. (2008) derive a putative orbital period for the wide 
binary of ~ 17 yr, implying an orbital velocity of approximately 
10- 15 kms"'. They also show that a circumbinary dust disk in 
such a system can explain the long-duration eclipse. A rotating, 
circumbinary disk would show an innermost velocity a little less 
than the binary orbital velocity, at some lOAU, in agreement 
with [O III] . The [O i] velocity corresponds to radii of a few hun- 
dred AU. 

The extreme Ha wings show the presence of a fast outflow, 
which in the presence of a stable disk will have a bipolar flow 
pattern. The dense wind shows the presence of outflow; we in- 
terpret this as a disk wind, caused by evaporation of the photo- 
ionized gas. The limit of stability of the gas disk is set when the 
thermal velocities exceed the rotational velocity. (Beyond about 
100 AU, only neutral gas or dust can be stable.) Evaporating gas 
can be accelerated by t he fast Ha wind, as a mass-loaded wind 
(iBorkowski et al.iri995h . The [Om] core emission shows a dou- 
ble peak, and this may indicate it forms in part in such a bipolar 
wind. 

A likely interpretation of the inne r nebula is, therefo re, in 
terms of an opaque disk seen edge-on (iHaiduk et aUlIOOSl) with 
ionized bipolar outflows. The infrared dust emission and neu- 
tral [Oi] line originate at the outer static boundary. The dense 
wind originates in the evaporating disk. At larger distances such 
an evaporating-d isk wind is expec ted to be essentially spheri- 
cally symmetric (lAIexandeijl2008h . The disk can be located in 
the same plane as the outer partial ring, at an inclination of 20 
degrees from the line of sight, with the bipolar flow directed to- 
wards the openings between the two arcs. 

6. Discussion 

6. 1 . The nebula 

The one-side structure of the wing in M2-29 is unusual 
but not uni que. It is duplicated in He 2-428, and in A 79 
(Rodriguez et al. 200 1[). and possibl v also in the symbiotic star 
V417 Cen dVan Winckel et alJll994 . In afl four cases, the struc- 
ture is essentially identical, allowing for the difference in view- 
ing angle. But no further object with a similar m orphology could 
be ide ntified in the lAC morphological catalog (iManchado et al.l 
|1996} ), confirming the rarity of this structure. 

Outer asymmetries can be formed by interaction with the in- 
terstellar medium, a nd these are a very common phenomenon 
dWareing et alj|2007h . In M 2-29, there is indeed a possible asso- 
ciation with interstellar gas at a similar velocity, and a possible 
ISM sweep-up should be explored. However, the good alignment 
of the outer ring with the equatorial plane of the nebula, the pres- 
ence of expansion, and the existence of a few very similar ob- 
jects, all suggest that the partial ring was part of the planetary 
nebula ejection. 

The arcs show the presence of a more (if not entirely) spher- 
ically symmetric component in the nebula, located within the 
outer ring. This suggests a scenario where the mass ejection was 
initially highly asymmetric, but became less so over time. The 
arcs could trace a temporary mass-loss enhancement, or they 
could be caused by the pressure enhancement from the onset of 
ionization. The compact core is evidence that a fraction of the 
ejected mass was retained close to the star 

The ejection of the partial ring can be related to an interac- 
tion with a companion with an elliptical orbit. Such an interac- 
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tion is discussed bv lSoker & Rappaportl ( l2001h . their scenario 4. 
If the orbit is such that interaction occurs at periastron only, am- 
plifying or initiating the mass loss at this point, the mass loss will 
occur preferentially in the direction of the velocity vector of the 
mass-losing star at periastron, and not at the systemic velocity. 
As the stellar radius increases, the binary orbit may circularize 
or the stellar mass loss may become self-sustaining: the mass 
loss will now occur everywhere along the orbit, centred on the 
systemic velocity. 

The interaction betwee n these two winds offset in velocity 
dSoker & Rappaporil2001h provides a natural explanation why 
an identical structure is seen in several objects, and seems the 
most plausible explanation for the outer wing, and the fact that 
the asymmetry occurred only during early mass loss. There is 
some evidence that V417 Cen, which shows a simil ar morphol- 
ogy, has an elliptical orbit dVan Winckel et al ] [T99l . 

6.2. Classifying tlie star 

M 2-29 is one of only three PNe known to sh ow dust eclipses (the 
others are: NGC2346 and CPD-56°803 2 jHaiduket al.ll200l 
ICohen et al.]l2002t iDe Marco etani2002h '). and one of two with 
a bin ary period of order 1 month (' NGC2346 has period of 16 
days (iHaiduk et al.ll2008tlde Marcolllool '). 

It is also one of few showing a dense, unresolved nebular 
core inside an extended FN. Other cases include EGB 6, a very 
faint and nearly circular PN around a binary composed of a very 
hot white dwarf and a visual point-lik e companion with dense 
emission-line nebula dBond et al.lll993h . jRodriguez et al.ll2001h 
finds a similar core in He 2 -428. BTel stars, a m ixed class which 
includes some young PNe ( lLamersetal.il 19981) . also show com- 
pact cores and Ha line wings similar to M2-29 ( IZickgrafl2003l) . 
But M2-29 lacks the low-ionization iron lines found in B[e] 
stars. A very late thermal pulse leads to the ejection of a new 
core (Ivan Hoof et al.ll2007l) . but this core is hydrogen-poor and 
different from what is seen in M 2-29. 

A relation to the dusty symbiotic stars has been proposed 
(iHaiduket al.i,200 8: Miszalski et al. 2Q09), based on the orbital 
period of M2-29. Symbiotic systems contain a mass-losing gi- 
ant, and a hot white dwarf which irradiates the wind: the ion- 
ized gas shows dense cores (e.g. M2-9, OH 23 1.8+04. 2). But 
although possibly related, M 2-29 does not appear to be a sym- 
biotic star It shows neither the Raman lines at 6830 & 7082 A, 
nor any indications for a cool (giant or yellow) companion in the 
spectrum. In the case of M 2-29, although there is evidence for 
two companions (Hajduket al. 2008), neither appears to be the 
required mass-losing giant. 

However, M 2-29 will likely be a symbiotic star at some time. 
If the companion is a white dwarf, the object would have been 
a symbiotic star or symbiotic Mira in the recent past, while los- 
ing mass. If the companion is main sequence, it will become a 
symbiotic once the companion ascends the giant branches (e.g. 
Ide Marcdl2009l) . 



6.3. Evolution of ttie disk 

Compact dust disks are common around post-AGB stars 
(Ide Ruvter et al.l l2006l) . and are increasingly being identi- 
fied around the cent ral stars of planetary nebulae: examples 



Table 3. Ages and dust masses of the disks in planetary nebulae 



Nebula 


Nebular Age 


Dust Mass 


Reference 




[yr] 


[Mq] 




CPD-56°8032 


102 


3 X 10-'' 


(2), (3) 


Ant nebula (Mz 3) 


10^ 


1 X 10"^ 


(1), (4) 


M2-29 


5x 10^ 


10-* 


This paper 


Helix nebula 


1.1 X 10" 


4 X 10^' 


(5), (6) 



References. (l)^hesneauetaL QOOd); (2)^hesneau^taL QOOT) 
(3) De Marco et al. (1997) : (4) Guerrero et al. (2004) : 
(5) Meaburn et al. (2008): (6) Su et al. (2007). 



The origin and evolution of these disks is not yet understood. 
They are interpreted as circumbinary disks which formed dur- 
ing the mass loss phase, or as debris disks either left over from 
the main sequence phase or built recently fro m tidally disrup ted 
smaller orbiting bodies as discussed in e.g. ISu et al.l (l2007h or 
iBrinkworth et al.l (l2009l) . The plethora of different classes of ob- 
jects with such disks (B[e] stars, symbiotic stars, RVTau stars) 
shows that they form relatively easily. Binary companions ap- 
pear required in all cases, but the orbital periods range from 200 
days (post-AGB binaries with disks) to 100 yr (such as M 2-9). 

Interestingly, similar dust disks have not been discov- 
ered around the post- common-envelope bin aries seen in 
about 15-20% of PNe dMiszalski et al l l2009l). with periods 
of hours to days. Although models ( Sandquist et al.l 119981 : 



are the ant ne bula (IChesneau et al 20071). the heli x nebula 



('Su et al .ll2007l). CPD-56°803 2 dCohen et all l2002l). M2-48 
(Philli ps & Ramos-Lariosl l2008l) and NGC 2346 (ICostero et al.l 
[19861) .' 



iNordhaus & BlackmanI l2006l) indicate that the common enve- 
lope is not necessarily completely ejected and there remains 
some gas around the system, there is no evidence that this resid- 
ual leads to a dust disk. The known dust disks are seen around 
longer-period systems which avoided a common envelope phase. 

M 2-29 is argued to trace the phase where the gas disk is lost 
and the dust disk remains. Rotation velocities far exceed thermal 
velocities while the gas is molecular or neutral. The dust wil l 
settle towards the mid plane at this time dPominik et al.ll2003h . 
But once ionization starts, the thermal velocities are comparable 
to the rotation velocities, and a disk wind initiates. In M 2-29, the 
disk wind began ~ lO-' yr ago, but has been decreasing over time, 
likely reflecting the decreasing gas reservoir. At the current time, 
some lO"'* Mq of gas remains. The current disk-wind mass-loss 
rate of 10"^ Mq yr"' gives a decay time of the disk similar to the 
age of the nebula. 

The future evolution is of interest to the disks seen around 
more evolved objects. The Spitzer spectrum with the strong sil- 
icate emission and PAH emission indicates that the dust is rel- 
atively hot. The dust will cool sharply once the star enters the 
white dwarf cooling track where it rapidly fades by a factor of 
100. The central star of the helix is such a cooling track star 
The cooler dust will emit mainly at longer wavelengths: the dust 
disk o f the helix is detected mainly at 24 jjm and 70 /im dSu et al.l 
l2007l) . Such cool disks can have gone largely unnoticed. 

Table |3] lists the ages and disk dust masses for four plane- 
tary nebulae where these have been measured. The ages are de- 
fined from the end of the AGB, i.e. the ejection of the main neb- 
ula. The dust masses are taken from the given references, where 
we assumed a gas-to-dust ratio of 100 for M2-29. The age of 
CPD -56° 8032 has likely been underestimated as the object has 
been known for longer than this dGilil900h . 

The Table indicates that the disk masses decrease sharply as 
the star ages. Although this should be taken with care because 
of the small sample and the possibly disparate origins of the ob- 
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jects, it appears possible that much of the dust mass is lost dur- 
ing the PN phase. Dust can be carried away in the wind from the 
evaporating gas disk, and unshielded dust may be des troyed by 
the intense UV radiation field. Chesneau et al. (I2006h suggests 
that an efficient dissipating process is occurring in the disk of 
CPD -56°8032. 

The original masses of the disks are unknown. However, 
disks around post-AGB stars show dust masses similar to 
CPD-56°8032: RU C en and AC Her sho w 5 x lO '* and 2 x 
10 "^Mq respectively ( iGielen et al.l l2007l) . These objects have 
stars which are still too cool to cause ionization, and have 
avoided ionization-driven dissipation. 

Thus, we suggest that the disk of the helix nebula may be 
the remnant of a post-AGB du st disk, rather than a debris disk 
as preferred by ISu et al 1 (I2007h . Once the star is on the cooling 
track, the disk becomes largely stable and the further dissipation 
time is set by Poynton-Robertson drag, with a time scale of or- 
der 10^ yr Thus, the disks may remain detectable for the order 
of 10*^ yr, long after any sign of the planetary nebula has disap- 
peared. 

Several disks have been discovered around white dwarfs (e.g. 
iGansicke et al.ll2006h . Surveys for such disks around the old- 
est white dwarfs have been unsuccessful (iKihc et al.ll2009h . but 
around 15% of local white dwarfs show metal-rich material in 
their phot osphere, indicativ e of accretion from a residual dust 
reservoir dSion et al ] l2009h . Although these may have a vari- 
ety of origins, they may include the progeny of the M 2-29-type 
disks. 



7. Conclusions 

We have presented a detailed study of one Bulge planetary neb- 
ula, with the aim to find evidence for disk evolution. Our main 
findings are: 

1 . M 2-29 has a likely distance consistent with membership of 
the Galactic Bulge. It shows absorption features of a cloud 
with high velocity, which we tentatively identify with the 
nearby HVC 002.9-06.2-088. 

2. The nebula is surrounded by a partial ring, slowly expanding 
at 12kms"'. Such a structure may be derived from interac- 
tion with the interstellar medium, but based on the expan- 
sion, and the alignment with the inner and outer nebula, we 
associate it with an early mass-loss phase from the central 
star A model where mass loss is triggered initially only dur- 
ing periastron of a companion in an elliptical orbit can ex- 
plain the structure: in this model, the ring was ejected with 
a preferential velocity in the orbital plane, reflecting orbital 
motion of the central star. 

3. We find evidence for on-going mass loss from the central 
source, decreasing with time, with a slow expansion veloc- 
ity of lOkms"', and with a current rate of ~ 10^^ Mqyt^^ . 
There is however no evidence for any mass-losing star which 
could be the source of this material. 

4. There is evidence for a very fast wind (~ lO^kms"') from 
the central star, based on the Ho- line wings. The [Oiii] also 
shows wings, but a factor of ten naiTower. 

5. The unresolved core with diameter < 250 AU shows strong 
forbidden and permitted emission lines with a range of ion- 
ization states, down to [Oi]. The velocity HWHM of the 
4363 [Oiii] line is 18kms"'; [Oi] is much narrower. Because 
of the very short dynamical time scale, and the velocity gra- 
dient with excitation, we interpret this as a stable, rotating 



system. The mass of the core is estimated as Mc ~ 10 Mq. 
The dust emission is located within this core. 

6. We interpret the on-going mass loss as a disk wind, driven 
by the ionization and heating of the gas. The thermal gas 
velocities are similar to the rotation velocity, making the disk 
susceptible to evaporation. The fast wind from the star may 
be accelerating gas shed from the disk, as seen in the [O iii] 
line wings. The mass and current mass-loss rate of the disk 
indicate an decay time similar to the age of the nebula. 

7. By comparing masses of dust disks in planetary nebulae of 
different ages, we show that the dust disks dramatically re- 
duce in mass as the nebula evolves, reflecting the disk evap- 
oration. The dust mass decrease from ~ 5 x 10"^ at the 
start of the PN phase, to ~ 4 x 10"^ Mq for the heUx neb- 
ula on the white dwarf cooling track where the luminosity of 
the star quickly drops by a factor of 100, reducing the rate at 
which the dust disk would be lost. 
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